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Introduction 
Mission Desi 


■ The interplanetary design problem is composed of both discrete and 

real-valued decision parameters: 

- Choice of destination(s), number of planetary flybys, identities of flyby 
planets 

- Launch date, flight time(s), epochs of maneuvers, maneuver 
magnitudes and directions, flyby altitudes, etc. 

■ For example, for a near-Earth asteroid mission, the designer must 

choose: 

- The optimal asteroid from a set of scientifically interesting bodies 
provided by the customer 

- Whether or not to perform planetary flybys on the way to the main belt 
and, if so, at which planets 

- Optimal trajectory from the Earth to the chosen asteroid by way of the 
chosen flyby planets 
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Automated M 
Optimal Cont 


■ Break the mission design problem into two stages, or “loops” 

- “outer-loop” picks sets of destinations, planetary flybys, sizes the 
power system, can pick propulsion system - a discrete optimization 
problem 

- “inner-loop” finds the optimal trajectory for a given candidate outer- 
loop solution - a real-valued optimization problem 

- For the outer-loop to work, the inner-loop must function autonomously 
(i.e. no human interaction) 
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Multi-Object 


■ The customer (scientist or project manager) most often does not want 
just one point solution to the mission design problem 

■ Instead, an exploration of a multi-objective trade space is required 

■ For a typical main-belt asteroid mission the customer might wish to see 
the trade-space of: 

- Launch date vs 

- Flight time vs 

- Deliverable mass 

- While varying the destination asteroid, planetary flybys, launch year, 
etc. 

■ To address this question we use a multi-objective discrete outer-loop 
which defines many single objective real-valued inner-loop problems 
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Outer-Loop T 



■ The outer-loop finds the non-dominated trade surface between any set 
of objective functions chosen by the user 

■ Non-dominated surface means “no point on the surface is superior to 
any other point on the surface in all of the objective functions” 

■ The outer-loop solver may choose from a menu of options for each 
decision variable 

■ The choices made by the outer-loop solver are used to define trajectory 
optimization problems to be solved by the inner-loop 
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• Break mission into a set of “journeys,” 
each of which in turn is broken into 
“phases” 

• The endpoints of a journey are chosen in 
the problem assumptions 


• The endpoints of a phase (i.e. a flyby 
target) may be chosen by the user or an 
Outer-Loop solver 
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Outer-Loop 


Launch Year 
Code Year 

0 2020 

1 2021 

2 2022 

3 2023 

4 2024 

6 2025 


Flight Time Upper 
Bound 


Code 

0 

1 

2 

3 

4 


# Years 

5 

6 

7 

8 
9 


First Asteroid 

ode 

Body 

0 

Ceres 

1 

Pallas 

2 

Juno 

3 

Vesta 

4 

Astraea 

5 

Hebe 


Second Asteroid 
Code Body 

0 Ceres 

1 Pallas 

2 Juno 

3 Vesta 

4 Astraea 

5 Hebe 


First Journey First Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


7 2026 

8 2027 

9 2028 

10 2029 


Second Journey Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


5 10 

7 11 

8 12 


6 Iris 

7 Flora 
(475 

^ choices) 


6 

7 


Iris 

Flora 

(475 

choices) 


First Journey Second Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 




Sample Mission 




Flight Time Upper 

Asteroid 1 

Potential Planetary 

Asteroid 2 

Potential 


Bound 


Flyby 1 


Planetary Flyby 2 

Code 

4 

0 

1 

1 

1 

Translation 

By 

Ceres 

Mars 

Pallas 

none 
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Multi-Objecth 


■ The outer-loop optimization problem is solved using a discrete multi- 
objective solver, in this case Non-Dominated Sorting Genetic Algorithm 
II (NSGA-II) 

■ NSGA-II finds the non-dominated front, surface, or hyper-surface 
between any number of objectives chosen by the user 


Initial generation 



Population evolves via 
genetic operators 

> 


Mass 



• • • 





TOF 


TOF 
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Inner-Loop M 


zation 


■ The inner-loop solves a real-valued trajectory optimization problem 
which is defined by each candidate solution to the outer-loop problem 

■ The inner-loop must function autonomously because the problems are 
generated in real time and there is no opportunity for human 
intervention 

■ The outer-loop is only as good as the solutions to the inner-loop 
problem, so the inner-loop must be robust 

■ A given run of the outer-loop may require hundreds or even thousands 
of runs of the inner-loop, so the inner-loop must be fast 

■ If the individual inner-loop runs are independent then many of them can 
be run in parallel 
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Multiple Gravity 
Maneuver (MGA 



3ace 
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Inner-Loop O 
Delivered Ma: 


■ Traditionally the objective function for a chemical mission is to minimize total Av because it is linear 

■ Av is used as an analog for mass via the exponential form of Tsiolkovsky’s rocket equation 

irif — m 0 e~ Av ^ c (1) 

■ However m 0 is actually a function of hyperbolic excess velocity C 3 , so just optimizing Av does not 
optimize spacecraft mass 

■ On the other hand, the derivative of (1 ) with respect to Av is very steep and therefore (1 ) is difficult for a 
gradient-based optimizer to handle. Instead we find a transformation of (1) works well: 

J — loglO(rrif ) = log 10 (m 0 e~ Av ^ c ) (2) 

■ Launch vehicles are modeled using a polynomial fit 

m 0 = (1 - a LV ) (i a LV Ci + b LV C$ + c LF Cf + d LV Cl + e LV C3 + / LF ) 
where o LV is a user-defined launch vehicle margin, zero for this presentation 
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Inner-Loop S 
Nonlinear Pr< 


Minimize f(x) 

Subject to: 
x lb <x< x ub 
c(X) < 0 

Ax < 0 

where: 

x ib , x ub are lower and upper bounds on the decision variables 
c(x) is a vector of nonlinear constraints 
Ax is a vector of linear constraints 

■ There are several third party solvers that do this (SNOPT, IPOPT, fmincon, 
vf 1 3AD) 

■ But all of these methods require an initial guess... 
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Inner-Loop 

Monotonic 



Leary, 2000 

Vasile, Minisci, and Locatelli, 2009 
Yam, di Lorenzo, and Izzo, 2011 
Englander (dissertation), 2013 
Casioli et al., 2013 
Englander and Englander, 201 4 


Improved from standard MBH by: 

1 . “Feasible point finder” aggregate penalty method 

2. Non-uniform (Pareto) perturbation step 

3. “Time-hop” operator (Casioli et al.) 
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Example: A 


Description 

Value 

Launch year 

outer-loop chooses in [2020, 20291 

Flight time 

outer-loop chooses in [3, 101 years 

Launch vehicle 

Atlas V 551 

Spacecraft Isp 

320 s 

Arrival condition 

insert into orbit at Jupiter 


a = 140RJ 


e = 0:91 

Number of flybys allowed 

up to 5 

Flyby targets considered 

Venus, Earth, Mars 

Outer-loop objective 
functions 

launch year 
flight time 
delivered mass 

Outer-loop population size 

256 

Outer-loop mutation rate 

0.3 

Inner-loop MBH run-time 

10 minutes 

Inner-loop MBH Pareto a 

1.3 
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Delivered mass 


Jupiter Mission: ■ 
First Generation ] 
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Delivered mass 


Jupiter Missic 
Final Generat 
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Jupiter Missio 
Example Traje 




8-year mission launching in 2021 5-year mission launching in 2026 
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y (km) 



Example: W 


Description 

Value 

Launch year 

outer-loop chooses in [2020, 20291 

Flight time 

outer-loop chooses in [3, 12] years 

Launch vehicle 

outer-loop chooses Atlas V 401 , 41 1 , 421 , 431 , 541 , 
or 551 

Spacecraft Isp 

320 s 

Penetrator mass 

20 kg 

Arrival conditions 


(first Journey) 

intercept with v M in [5.0, 10.0] km/s, 0m liminatinn ^ 70° 

(second Journey) 

rendezvous 

Number of flybys allowed 

2 in each Journey 

Flyby targets considered 

Venus, Earth, Mars 

Outer-loop objective 
functions 

launch year 
flight time 
delivered mass 
launch vehicle choice 

Outer-loop population size 

256 

Outer-loop mutation rate 

0.3 

Inner-loop MBH run-time 

10 minutes 

Inner-loop MBH Pareto a 

1.3 
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Delivered mass to final target (kg) 


Whack-a-RocI 
First Generati 
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Launch vehicle preference 


Delivered mass 


Whack-a-Rock 
Final Generate 
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Launch vehicle preference 



Whack-a-Roc 
Example Traj< 


Event # 2: 
chemical burn 
deep-space 
4/30/2029 
Av = 0.723 km/s 
m = 3157 kg 


Event # 4: 
chemical burn 
deep-space 
3/13/2033 
Av = 0.224 km/s 
m = 2920 kg 


Event # 5: 
unpowered flyby 
Earth 
12/3/2033 
v x = 4.560 km/s 
DEC = 64.9 
altitude = 20409 km 
m = 2920 kg 


Event # 3: 
intercept 
162173J.999JU3 
1/10/2030 

v^ = 8.122 km/s 

DEC = 47.1° 
m = 3157 kg 



Event # 7: 
rendezvous 
162173JL999JU3 
1/23/2037 
Av = 0.018 km/s 
m = 2859 kg 


Event # 1: 

launch 

Earth 

10/27/2025 
C3 = 8.026 km 2 /s 2 
DLA = 22.4 
m = 3975 kg 


Event # 6: 
chemical burn 
deep-space 
6/1/2036 
Av = 0.048 km/s 
m = 2876 kg 



Atlas V 421 , 11 .25 year flight time 


Atlas V 551 , 2.45 year flight time 


NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


21 * 


y (km) 


Conclusions 




■ The chemical interplanetary mission design problem may be posed as a multi -objective hybrid 
optimal control problem 

■ The combination of a multi-objective discrete NSGA-II outer-loop with a MBH+NLP inner-loop is 
a very powerful way to explore a mission trade space in an efficient, automated manner 

■ The algorithm described here is a valuable force-multiplier for interplanetary trajectory design 

We can now study multiple mission design cases simultaneously, limited only by available 
computing power 

Mission design engineers can now spend more time with the customer and with spacecraft 
hardware engineers so that we can fully understand the scientific and engineering context of 
our work 

Good mission ideas are much less likely to be rejected due to lack of time to work on mission 
design, and bad ideas are much more likely to be rejected before they consume too many 
resources 

■ Skilled analysts are expensive. With a multi-objective HOCP automaton, analysts can focus on 
understanding the customer’s needs and the spacecraft’s capabilities and also detailed design 
work, leaving repetitive tasks to the computer 
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EMTG is available open-source at 
https://sourceforae.net/proiects/emta/ 
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